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Abstract. Let C be an irreducible projective curve of degree d in P n (K), 
where K is an algebraically closed field, and let / be the associated homoge- 
neous prime ideal. We wish to compute generators for /, assuming we are given 
sufficiently many points on the curve C. fn particular if / can be generated 
by polynomials of degree at most m and we are given md + 1 points on C, 
then we can find a set of generators for /. We will show that a minimal set 
of generators of / can be constructed in polynomial time. Our constructions 
are completely independent of any notion of term ordering; this allows us the 
maximal freedom in performing our constructions in order to improve the nu- 
merical stability. We also summarize some classical results on bounds for the 
degrees of the generators of our ideal in terms of the degree and genus of the 
curve. 



1. Introduction 

Let C be an irreducible projective curve of degree d in P™(K), where K is an 
algebraically closed field, and let / = 1(C) be the associated homogeneous prime 
ideal of V = K[x , ■ • ■ > x n ] consisting of all the polynomials vanishing on C. We 
wish to compute generators for /, assuming we are given sufficiently many points on 
the curve C. In particular if / can be generated by polynomials of degree at most 
m and we are given at least md+ 1 points on C, then we can find a set of generators 
for /. It is a simple consequence of Bezout's theorem that any polynomial of degree 
k which vanishes on more than kd points of C must be contained in /. Although 
the number of monomials in n + 1 variables of degree at most m is not polynomial 
in both n and m, we will present a process which constructs generators degree 
by degree and results in a polynomial time algorithm for computing generators 
for /. Polynomial time algorithms for computing Grobner bases of ideals of affinc 
points were presented in MB , and then extended to minimal generators of ideals 
of projective points in [MMM . These algorithms require exact arithmetic, and 
assume a term ordering is given. We present new algorithms which are completely 
independent of any notion of term ordering. We believe that this flexibility is 
necessary when working with approximate coefficients. 
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Given a homogeneous ideal /, there are many different choices for monomials 
representing cosets of V/I, i.e. for a complement of /. It is well known that the 
natural coset representatives associated with Grobner bases do not remain stable 
with respect to small coefficient perturbations of the ideal generators. Border bases 
were introduced to help overcome this problem (jKRj). Given a fixed choice of 
complement for a zero-dimensional ideal, its border basis is uniquely determined in 
contrast with Grobner bases where the complement is uniquely determined by the 
given term ordering. Border bases are usually defined for zero-dimensional ideals, 
which guarantees a finite basis. We extend the definition to homogeneous ideals 
of any dimension, but bound the degree in order to preserve finitcness. In much 
of the literature on border bases, the complements are required to be closed under 
division by variables. This makes complements of border bases more similar to 
complements of Grobner bases which also have this property. In particular this 
is done in [HKPPj . therefore their algorithms need to explicitly decide whether 
or not candidate leading monomials have coefficients which are so small that they 
should be treated as zero. As suggested by Mourrain and Trebuchet ( |MT| ) . we 
only require the complement to be connected to 1. This means we require that each 
complement monomial of degree i is a multiple of some complement monomial of 
degree i — 1. This extra flexibility in the choice of complement monomials means 
that we can use standard numerical software like the QR algorithm with column 
pivoting (QRP) ( |GVLj ) to choose our complement in each degree. One could 
define a complement to be any set of representatives for V/I, but with this definition 
we would not be able to obtain algorithms which are polynomial in both the degree 
of the curve and the number of variables. In particular, requiring the complement 
to be connected to 1 implies a strong condition on the syzygy module. We show 
that the syzygy module for a vector space basis of a homogeneous ideal whose 
complement is connected to 1 is generated by vectors whose entries have degree 
at most one, generalizing the result of Mourrain and Trebuchet for border bases 
of zero-dimensional ideals. These special generators of the syzygy module can be 
used to obtain a polynomial time algorithm for constructing minimal generators 
for our ideals. The algorithms developed by Cioffi [C\ have a similar complexity in 
the case of exact coefficients but her use of Grobner bases requires a term ordering 
which determines a unique complement which may be numerically unstable when 
working with approximate points. 

The motivation for this paper came from the desire to be able to compute the 
generators of space curves starting from numerical software which generates points 
on curves, in particular computing points on canonical or bicanonical models for 
Riemann surfaces presented as Fuchsian groups ( |GSSTj ). Our intended appli- 
cations differ from that of [HKPPj and |AFT] since we assume that the points 
defining our curve are generated by a numerical algorithm whose errors tend to 
be very small as opposed to empirical measurements whose errors could be much 
larger. Thus our goal is to present algorithms based on numerically stable con- 
structions like SVD for computing ideal generators and QRP for deciding which 
monomials represent the complement of our ideal. 

In section 2 of this paper we derive some general properties of border bases and 
complements for general homogeneous ideals. Using these properties we present 
a polynomial time algorithm for finding a minimal basis for a homogeneous ideal. 
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In section 3 we specialize to ideals of curves given by points, and use point eval- 
uation matrices to complete the task of computing border bases for homogeneous 
ideals of curves. Assuming exact arithmetic of unit cost, we also provide an overall 
complexity analysis of the border basis algorithm and the minimal basis algorithm. 
In section 4 we consider the situation of approximate points and show how our 
algorithms can be adapted to use standard numerical software, where we allow nu- 
merical algorithms like the QRP to choose the complement monomials in order to 
help improve the numerical stability. In general the stability also strongly depends 
on how the points are distributed, and we feel this is an interesting problem for 
future research, along with the possibility of using other approaches such as interval 
arithmetic. 

In order to use Bezout's theorem, we assume we have at least a bound on the 
degree of our curve. In the last section we summarize some classical results on 
bounds for the degrees of generators of our ideal in terms of the degree and genus 
of the curve. In particular, as shown by Petri ([P]), a non-hyperelliptic canonical 
curve of genus g > 4 can be generated in degrees 2 and 3. We also recall the 
completely general result of Gruson-Lazarsfeld-Peskine f |GLPj ) which shows that 
any non-degenerate curve of degree d in P™(K) can be generated in degree d — n + 2. 

Sometimes we have additional information about the nature of the curve whose 
points we are given. For instance, if we know the Hilbert function, the rank of 
our point evaluation matrices is explicitly given instead of being determined by 
examining its singular value spectrum. 

2. Border bases for homogeneous ideals 

Let K be an algebraically closed field. For any s £ N, let T s be the set of all 
terms of degree s in V = K[xq, . . . , x n ] and let V s be the vector subspace of V 
generated by Ts- Recall that dim^ s = ("+ s ). We will also set V< s = ®i =0 Vi. 

Notation 2.1. We will use the following notation: 

(1) For any subset Y C P n (K) denote by T{Y) the radical homogeneous ideal 
of V consisting of all the polynomials vanishing on Y . 

(2) For any homogeneous ideal I in V , let I a = I D V s (so that / = ® s >oI s ) 
and let I< a = I C\V< S . 

(3) For any ScP, denote by I{S) the ideal generated by S. 

(4) For any S C V s , denote by (S) the vector subspace of V s generated by S. 

(5) For any S C V s , let S+ = U"=oOj S) C V s +x- 

(6) If a = {a x , ...,a h )£ K h and T = [Fx,..., F h ] is a list of polynomials, we 
set 

a • T = axFx + . . . + a h F h . 

(7) For any finite set A, we denote by \A\ its cardinality. 

Definition 2.2. Let J be a proper homogeneous ideal in V and s £ N. Let 
A/o = {1} and, for each k — 1, . . . , s, assume that A/& is a set of monomials in Tk 
such that 

M k C A^lj and V k = J k ®{N k ). 
We call J\f — {A/o, • ■ • , AA S } a complement of the ideal J up to degree s. 

Remark 2.3. Let J\f = {A/o, ■ ■ • ,A/" S } be a complement of a proper homoge- 
neous ideal J up to degree s. Then: 
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(1) the condition Afk C Af^_ x implies that Af is connected to 1, i.e. for each 
to G Af there exist variables Xi t , . . . , Xi k such that m ~ Xi x ■ . . . ■ Xi k and 
Xi x ■ . . . ■ Xi, S Af for each j < k, 

(2) from the definition it follows that, for k = 1, . . . , s, 

(Afcj) = (J k n (A/"+_ x » © (A4). 

□ 

Lemma 2.4. Let J be a proper homogeneous ideal ofP. Assume that Afk-i <!= 
7fe-i and TVfc C 7fe are sets o/ monomials such that 

(a) TVi = Jfe-i © (M;-i> 

(b) Wt-i) = V* n(AT+_ 1 )) ©(A/"*). 

TTien 

(1) J fc = <j+ ^ + (J fe n {Af+_ x )) 

(2) = J fe ®(A4). 

Proof. By hypothesis (a), we have 

n = (H-i) = (4-x) + (K-x) 
and, since Jfe 2 (J^_i), we have also 

^ = p fc n Jfe = (Jfeli) + (Jfe n M+J), 

which proves (1). Hence, by the previous relations and hypothesis (b), we have 

n = <J fc ti> + Wti) = + ( J fc n + M> = J * + Wt>- 

On the other hand, again by hypothesis (b) we have that (Afc) C (Aft-.) and hence 

Jfe n (A/3k> c Jfe n (A/"+_ x ) n (Af k ) = {0}, 

which completes the proof of (2). □ 

Remark 2.5. It is always possible to choose a complement up to any fixed 
degree for any proper homogeneous ideal J of V incrementally. Namely, if Af = 
{Afo, ■ ■ ■ , A/fc-i} is a complement of J up to degree k — 1, it is sufficient to choose 
a set TVfc C 7^ such that (Af^) = (J fc n (Af^)) © (A/fe): then Lemma \2A\ assures 
that Vk — J k © (Afc) and hence that Af = {Afo, ■ ■ ■ ,Afk} is a complement of the 
ideal J up to degree k. 

Moreover, again by Lemma l2.4[ once one has a complement Af of J up to any 
fixed degree s, one can get a set of generators of the ideal I{J\, ■ . ■ , J s ) as the union 
of sets of generators of Jfe fl (Afj^_i) for k = 1, . . . , s. 

Definition 2.6. Let J be a proper homogeneous ideal in V and assume that 
Af = {Mq, ...,Af s } is a complement of J up to degree s. 

(1) For all k = 1, . . . , s let (cW)fc = A/^ \ A4; the elements in (<9A/jfc will 
be called border monomials in degree k. 

(2) For each m £ (dAf)k let ^(m) be the unique polynomial in (A4) such that 
to + ip(m) £ Jfe; we will call the homogeneous polynomial m + if)(m) the 
border polynomial associated to m. 

(3) If Bk denotes the set of all border polynomials of degree k, the set B = 
B\ U . . . U B s is called the border basis of J up to degree s associated to Af. 
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Remark 2.7. This notion of bounded degree border basis applies to arbitrary 
homogeneous ideals and coincides with the classical notion of border basis ([KR ) 
in the case of homogeneous zero-dimensional ideals, provided we choose s to be 
larger than the maximal degree of any monomial in the (finite) complement J\f of 
J. □ 

Proposition 2.8. Assume that N is a complement of a proper homogeneous 
ideal J up to degree s and let B = B\ U . . . U B a be the associated border basis. 
Then 

(1) Bk is a basis of the vector space </& n (A/iti ) for each k = 1, . . . , s, 

(2) B is a set of generators of the ideal I(Ji, . . . , J s ). 

PROOF. (1) Let b G J k n (A/^). Since A/jL x = (cW) fc U Af k , there exist 
Oj, bi G K such that 

b= 22 aimi+ ^ hrrii- 

m;6(aV)fc HH6A/* 

For each G {dN)k let if) (mi) be the unique polynomial in (A4) such that m, + 
il)(rni) G Jfc. Then we can write 6 = it + v with 

u= 2J a,i(mi + ip(m,i)), v= 2J - a^mi). 

Note that u € Jfc D (Sfc) and i> G (A4)- Since 6 G Jfc, then we have that u = & — it 6 
^fc n (A/fe) = {0}. Thus fe = u and hence 6 G (Bk). 

(2) follows immediately from (1) and Remark 12.51 □ 

We now suggest a simple method to construct recursively both a complement 
and a border basis of a homogeneous ideal J up to any fixed degree. Even if 
the ideal may not have an explicit representation, if we assume the capability of 
computing a basis of its intersection with a vector subspace generated by a finite 
set N of monomials of the same degree, we will be able to compute a border basis 
for J up to any fixed degree. We will denote this condition by saying that the ideal 
is represented by the function ComputeBasisj , which, for any such N, returns 
a basis for the intersection J n (N). In the next section we will see that such a 
function can be easily computed for ideals of points. 

In the description of the algorithms we will use the following notations: 

a. If v is a polynomial and S = {ni, . . . ,n t } is a set of monomials, then 
cocffs(i>, S) will denote the vector (ai,...,a t ) of the coefficients of the 
monomials of S in v. 

b. Given a matrix A, we will denote RRE(A) — (E, E) where 

- E is the completely reduced row echelon form of A (i.e. each pivot 
is equal to 1, and in each of the columns containing a pivot all the 
elements different from the pivot are zero) 

- E is the set consisting of the indexes of the columns containing the 
pivots of E . 

Algorithm BorderBasisWithComplement 
Input: 

- a function ComputeBasisj representing a homogeneous ideal J 

- s GN 
Output: 
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- {A/o, • • • , A/" s } a complement of J up to degree s 

- {Bi, . . . , B s } the associated border basis. 

Procedure: 

- A/o = {1} 

- for fc = l..s repeat 

— construct the set of distinct monomials = {mi, . . . ,mt} 
J^k-i := { x i m I 1 < i < n,m £ A/fe_i} 

Vfc := ComputeBasisj{M^_ l ) 
q-=\Vk\ 

f-=WLi\ 

— compute the q x t matrix with rows coeffs(«, A/jL-,) for »e Vi 
A := matna;(coeffs(w, A/^Lj) | w e V&) 

(E 1 , S) := RRE(A) 

— the monomial with index not in E are put in A4 

A4 := {m,j | j £ £} 

— every row represents an element of Bk 
®k ■= {J2j e i,j m j I* = !>•••) ?} 

Proposition 2.9. Given a proper homogeneous ideal J in V represented by a 
function ComputeBasisj and s 6 N, the algorithm BorderBasisWithComplement 
constructs a complement and a border basis for J up to degree s. 

Proof. At each step, the rows of E correspond to polynomials of the form 
m + ip(m), with m A4 and ^(m) G A4, which are a basis of Ju n (M^Li)- 
Moreover , since the monomials in A4 correspond to the non-pivot positions, by 
construction we have that (Af^-i) = (Jk H (M^^) © (A4). 

Then by Lemma 12.41 we get that Vk = Jk ffi (A4); hence {A/"o, ■ • ■ , A/fc} is a 
complement of J up to degree k and £>i U . . . U Bk is the associated border basis up 
to degree k. □ 

By Proposition 12.81 a border basis B of J up to degree s is a set of generators 
of the ideal I(Ji, . ■ . , J s ), but in general it is not minimal. We will see how one can 
eliminate redundant polynomials in B so as to obtain a minimal set of generators 
of that ideal. 

Proposition 2.10. Let J be a proper homogeneous ideal in V. Assume that 
J\f = {A/o, ... ,J\f s } is a complement of J up to degree s and let B = Bi U . . . U B s 
be the associated border basis. Then: 

(1) The ideal L = 7(£>,A/" S + ) is homogeneous and zero-dimensional, and B U 
7V S + is the border basis of L associated to its complement TV. 

(2) The module Syz(B, W s + ) is generated by vectors whose entries have degree 
at most 1. 

Proof. (1) Since L s — J S: we have that *P S — L s © (A/^). Moreover, since 
A/j+ C L, choosing Af s+ i = we see that the triple L,J\f s ,J\f s+ i satisfies the hy- 
potheses of Lemma \2 .41 (we set (0) = {0}). Thus we get that V s +i — L s +i\ so the 
ideal L is zero-dimensional. If we set A/} = for all j > s + 1, then TV = {A/jjjgN 
is a complement of L. Moreover, since (<9A/")s+i = A/" s + , if we let B s+ i = A/" s + , then 
B U A/" s + is a border basis of L. 
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(2) Since Nk+i Q Nu for each k and Ao = {1}, we have that TV" is connected to 
1 (see Remark l2.3p . Thus it is possible to apply Theorem 4.3 in }MTj . which implies 
that the syzygies among the elements of a border basis of a zero-dimensional ideal 
with a complement connected to 1 can be generated by syzygies whose coefficients 
have degree at most 1. □ 

By changing bases in each constant degree subspace we obtain the following 
more general result: 

Proposition 2.11. Let J be a proper homogeneous ideal in V. Assume that 
J\f = { Ao , • ■ • , Af s } is a complement of J up to degree s . Let V = Vi U . . . U V s where 
Vfc is a basis of Jk f~l (NjJLi) f° r eacn ^ = 1, • • • , s. Then the module Syz(V,Af+) is 
generated by vectors whose entries have degree at most 1. 

The following corollary shows that the redundant elements in V s can be ex- 
pressed as a combination of elements in V s _i and the other elements in V s : 

Corollary 2.12. Under the hypotheses of Proposition 12. 11[ let / £ V s and 
denote W s = V s \ {/}. If / £ I(Vi, . . . , V._i, W.), then / e (V+_i, W.). 

Proof. By hypothesis there exists a syzygy among the elements of V = Vi U 
. . . U V s such that the coefficient of / is a non-zero constant. Hence by Proposition 
12.111 there exists a homogeneous generator of Syz(V, A/" s + ) whose entries have degree 
at most 1 and where the coefficient of / is a non-zero constant. Since deg / = k 
and the coefficient of / is constant, in this generating syzygy only the elements of 
Vfc_i U Vfc can have non-zero coefficients. □ 

We now describe two methods to construct a minimal set of generators of the 
ideal I(Ji, . . . , J s ) depending on whether we start with a set of generators which 
form a border basis or not. 

Proposition 2.13. Let seN and let J\f = {No, . . . ,N S } be a complement up 
to degree s of a proper homogeneous ideal J in V . Given {Vi, . . . , V s } where Vfc is 
a basis of J (~l (A/^Lj), then for each k = 1, . . . , s it is possible to construct a set of 
polynomials Qk C Vfc such that Gi U . . . U Q s is a minimal set of generators of the 
ideal I{J\, . . . , J s ). 

Proof. Let Q\ = Vi and assume that Q\ U . . . U Gk-i is a minimal set of 
generators of I(Ji, . . . , Jfc-i) with Qi C Vj for i = 1, . . . , k — 1. 

Note that a polynomial / £ Vfc is redundant w.r.t. Gi U . . . U Gk— l U Vfc if and 
only if it is redundant w.r.t. Vi U . . . U Vfc . Thus, by Corollary |2.12l it suffices to look 
for linear relations among the elements of V£_ 1 U Vfc and for a set Gk Q Vfc such that 
( v k-i) + ( v fe) = ( v fc-i) © (Ok)- Hence it is sufficient to find a basis of (V^_ x ) H (Vfc), 
extend it with elements wi, . ■ ■ ,w t to a basis of (Vfc) and define Gk = {tvi, ■ ■ ■ ,w t }- 

In order to compute the intersection (Vti) ^ (^fe) consider the monomial basis 
S^S 1 US 2 US 3 oi (y+_ t ) + (Vfc), where Si = [BJf)t-i & - (^AOfc and 

S3 = A4. Observe that AA+_ 1 = (9A0fc UA4 and V^_ x C ((dU)^) + (A/^A. Let 
Sj = |Sj|, for z = 1,2,3, and I = |V fc f _ 1 |; with this notation |Vfc| = \(dAf)k\ = s 2 . 

Let 17 be the matrix whose columns contain the coefficients of the polynomials 
of V^-i u ^fc with respect to S. Thus [/ is a (si + S2 + S3) x (Z + s 2 ) block matrix 
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of the form 



u= \ u 2 



u. 



u 5 J 



and, if we denote by tt 2 : K l x K S2 — > K s ' 2 the projection on the last s 2 coordinates, 
the vectors of n 2 (KerU) are the coordinates (w.r.t. Vk) of the vectors of (V A L 1 ) n 

(V fc >- 

In order to compute Ker U we can reduce ourselves to consider the matrix 

U 








u 2 


Us 



Namely, since S3 = Nu, if Uv = 0, then Uv G J& n (A4) = {0}, hence Ker?7 = 
Kerf/. 

In order to finish the construction it is then sufficient to reduce to echelon form 
the matrix whose rows are generators of 7r 2 (Ker U) : the indexes of the columns 
without pivots correspond to the elements in Vk to select for constructing Q k ■ □ 

The proof of the previous proposition guarantees the correctness of the follow- 
ing: 

Algorithm MinimalBasis 
Input: 

- s e N 

- {Wo, . . . ,ftf s } a complement of a homogeneous ideal J up to degree s 

- {Vi, . . . , V s } where Vk is a basis for J n (Nk-i) 
Output: . . . , Q s } where: 

Qk C Vk for each k G {1, . . . , s} 

- the polynomials in Qi U . . . U Q s are a minimal set of generators of the ideal 

I(Jl,.;Js) 

Procedure: 

-Qi=Vi 

- for k = 2..s repeat 

S 2 ■■= (dM)k 

i:=|V fc + _il 
si := \Si\ 
s 2 := \S 2 \ 

— construct the (si + s 2 ) x (I + s 2 ) matrix with columns the s\ + s 2 

— coefficients w.r.t. Si U S 2 of the polynomials v s V^_ l U Vk- 
JJ := matrix (coeffs(u, Si U S 2 ) v G U V k ) 

— compute the intersection (V fe " t l 1 ) fl (Vfc) 

K := Ker((7) 
dfc := \K\ 

— construct the dk x s 2 matrix with rows the last s 2 entries of the 

— vectors in K. 

MK := matrix(jT 2 (v) \v G K) 
(RMK, S) := RRE(MK) 

— the polynomials in Vk with index not in S are put in Qk 
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Gk := {vj G V fe \j £ £} 

In the case when we start with a border basis, we can improve our previous 
construction. The computation of the generators of the intersection (V r fc l l 1 ) n (Vfe) 
can then be accomplished with only some column subtractions, and the construc- 
tion of each level of the minimal basis can be completed with one column echelon 
reduction of an s 2 x (I — s\) matrix. The new algorithm is based on the following 
Proposition whose proof is an immediate consequence of the properties of a border 
basis. 

Proposition 2.14. Let s G N and let J be a proper homogeneous ideal in 
V . Assume that M — {A/"o, ...,Af 8 } 1S a complement of J up to degree s and let 
B = B\ U . . . U B s be the associated border basis. With the notation of the previous 



proof, if U 



U 2 







U 3 



we have: 



(i) after reordering the elements of B k we can assume that the s 2 x s 2 block 
U3 is the identity matrix J S2 

(ii) each element in U\ is either or 1; more precisely each row in U\ contains 
at least one element equal to 1 and each column in U\ contains at most 
one element equal to 1 

(iii) by means of finitely many subtractions performed on the left I columns of 
U we can reduce U to the form 

U = 



I S1 








Pi 


p 2 


Is, 



(iv) the set of the columns of the s 2 x (I— Si) matrix P 2 is a basis of n 2 (Ker U) = 
7r 2 (KcrC/) 

(v) U has full row-rank and so dimKer U — dimKer U = I — s\ 

(vi) reducing P 2 to echelon form by column operations, the pivots indicate the 
redundant elements in Bk- 



3. Curves given by points 

A natural application of the results of the previous section is the construction 
of the ideal of an irreducible projective curve starting from the knowledge of a finite 
set of points on it. Let us recall the following classic result: 

Proposition 3.1. Assume that C is an irreducible projective curve in P™(K) 
of degree d. Let 1Z — {i?i, . . . , Rh} be a set of points on C. 

(1) For all s G N such that h > sd, we have 1(C) < s = 1(U)< S . 

(2) If 1(C) can be generated by polynomials of degree at most m and h > md, 
thcnX(C) =V-X(U)< m 

PROOF. (1) It suffices to prove that l{K) k C 1(C) k for all k < s. If / G l(TZ) k , 
the polynomial / vanishes on h > sd > kd = degf ■ degC points. Since C is 
irreducible, by Bezout's Theorem the hypersurface V(f) contains C, i.e. / G l(C) k - 

(2) By hypothesis 1(C) = V -l(C)< m , thus the result follows immediately from 
(1). □ 

The previous result allows us to reduce the construction of 1(C) to the compu- 
tation of a set of generators for the ideal J = 1(H) where 1Z — {Ri, ■ . ■ , Rh} is a set 
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of h points in P" (K) . By Proposition 12.81 this can be done by computing a border 
basis of J. In the case of an ideal of points, we are able to compute Jk l~l (NjT^) 
using the point evaluation maps. 

Assume that we have computed A4-i and let Afj[_i = {mi, . . . , mt} C Vk- 
Consider the h x t evaluation matrix 



mi(.Ri) ... m t (Ri) 



m 1 (R h ) ... mt(Rh) 



where, if Ri = [fYo, ■ ■ ■ , ?"i,nL by mj{Ri) we mean mj(ri t a, ■ ■ ■ j ^i.n)- Note that the 
rank of the matrix Mr and its null-space Ker M-ji does not depend on the chosen 
representation of the points in the projective space. Note also that each vector in 
Ker Mr is the vector of the coordinates of a polynomial in Jk (1 (Nk-i) w - r -t- the 
basis {mi, . . . ,mt}. In particular dimKerMR = dim(Jfe n {NuLi))- 

Performing Gaussian elimination by rows, followed if necessary by a permuta- 
tion of the columns (which corresponds to a permutation of the basis {mi, . . . , mt}), 

In this way the columns of [ — = ) are 



Ir 


A ' 





, 



we can assume that Ai-jz = . , . • •>.., , 

V U U J \ J-t—r 

a basis of the null-space Ker Mr. If we choose as A4 the first r monomials in the 
permuted basis, we have the null-space in border form, which gives us Bk- 

We now want to estimate the complexity of our procedure to compute a minimal 
set of generators up to degree s of the ideal J of h distinct points. 

Our algorithm first computes a border basis up to degree s, then minimizes 
this basis removing redundant elements. The basic tool is Gaussian elimination; 
recall that the complexity of Gaussian elimination performed on a to x p matrix is 
0(mpmiri(m,p)). 

As for the first phase to compute the border basis, the fc-th step of the recur- 
sive procedure described above to compute Nk and Bk requires to perform Gauss- 
ian elimination on the h x t matrix Mr. As already observed, dim Ker Mr = 
dim(Jfc n (A/'jLj)); hence, by Remark 12.31 dimKer Mr = t — dim(A4}- In particu- 
lar dim(A4) = |A/fc| = rk M-ji < h. Thus for each i — 1, . . . , s we have that \Afi\ < h 
and hence \Af^~ \ < (n + l)h. Therefore the complexity of each step of the algorithm 
is 0(nh 3 ). As a consequence the complexity of the recursive algorithm in s steps 
to compute a border basis of I(Ji, . . . ,J S ) is 0(snh 3 ). 

As for the minimizing phase outlined in Proposition 12 . 141 note that for each k 
we have that \B k \ = \(dAf)k\ <(n + l)h, (B^J < (n + l) 2 h and hence 

\Bk\ = {B+^UBkl < (n+l)(n + 2)/i. 

Moreover the distinct monomials appearing in the polynomials of Bti U Bk form 
a subset of anc ^ therefore they are at most (n + l) 2 h. 

Since the number of the left columns in U is | < (n + l) 2 h, when we reduce 

the matrix U to the form U by means of subtractions among the left columns, we 
need at most (n + l) 2 h column subtractions. The length of each of these columns 
is at most (ra + l) 2 h, thus the complexity of the reduction of U to U is 0(n 4 h 2 ). 

The last part of the minimizing phase requires to reduce the matrix Pi to 
echelon form by column operations. The number of rows of Pi is equal to | {dN)k \ < 
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(n + l)ft, while the number of its columns is < (n + l) 2 h. Thus the complexity of 
the algorithm to reduce Pi is 0(n 4 h 3 ). 

Hence the complexity of the algorithm outlined in Proposition 12 . 141 to compute 
a minimal set of generators for J< s is 0(sn ft ). 



In this section we consider the situation where the points JZ = {Ri, ■ ■ ■ , Rh} 
on the irreducible projective curve C in P™(C) are given only approximately. In 
this case we can perform the computations needed for the described procedure by 
replacing Gaussian elimination by more suitable and numerically stable tools. 

As observed in Section [21 if one is only interested in computing a set of gener- 
ators of Z(Ji, . . . ,J S ) (not necessarily a border basis up to degree s), it is sufficient 
to compute a basis of Jfc Pi (•A/ti) and compute a complement Aff- for each k (see 
Remark |2.5|) . The former task corresponds to computing a basis of the null-space 
of the h x t matrix 



where {mi , m t } = A/^Lj. 

A numerically stable way to compute both the rank and an orthogonal basis of 
Ker M-ji is the SVD-algorithm which assures that one can find a unitary ft x ft matrix 
U, a unitary t X t matrix V and an ft x t real matrix £ such that Mn. = UHV ; the 
elements ay of the matrix S are zero whenever i ^ j and for i = 1, . . . ,1 = minjft, t} 
we have ai.i > . . . > on > 0. 

Either we know the rank of Mji (see for instance Proposition I5.8|) or we can 
examine the singular values ai of £ in order to obtain a rank determination as in 
|G VL] . In any case if rkMR = r, then by the properties of the SVD decomposition 
the last t — r columns of V are an orthogonal basis of Ker M-ji . 

In order to compute A4 and continue to the next step, we take advantage of 
the stability properties of the QRP-algorithm which, given a matrix M, constructs 
a unitary matrix Q, a permutation matrix P and an upper-triangular matrix R = 
( Ri | i?2 ) such that S = QRP. The permutation matrix P exchanges columns 
in order to improve the condition number of the matrix R\ . If M has full row-rank, 
then i?i is invertible; otherwise it is possible to use the diagonal elements of R\ to 
make a rank determination of M . 

In our case we apply the QRP-algorithm to the (t-r)xt matrix S whose rows 
are the last t — r columns of V; recall that the columns of 5* are indexed by 
In the decomposition S — QRP the columns of R are a permutation of the columns 
of S and the monomials corresponding to the columns of R\ will be chosen to be 
border monomials, while the monomials corresponding to the columns of i?2 will 
be chosen as the complement A4 . Observe that the rows of R correspond to a new 
basis for J k n (A/^ l l 1 ). 

If we want to compute a border basis of I{J\ 1 ■ ■ ■ , J s ), we can compute the 
matrix R^R = ( I | i?^ 1 i? 2 ) whose rows correspond to a basis £>& of KerM-^ 




consisting of border polynomials. 

Otherwise, if we want to compute a minimal set of generators of I(Ji, ■ ■ ■ , J s ), 
we can proceed as described in the algorithm MinimalBasis using SVD to compute 



4. Curves in P"(C) given by approximate points 
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kernels and QRP to select stable pivot columns. Using the notation of Proposition 
12.131 the first step is to compute a basis of Ker U. In order to do this, we use an 
SVD construction taking into account that, by Proposition 12.141 (v) dim Kerf/ = 
dim Ker U = I — s\. We then apply the QRP-algorithm to the matrix N whose rows 
contain the projection by 7r 2 of a set of generators of Ker U and whose columns 
are indexed by the generators of J k n (7V / ,jl 1 ). We thus obtain N — Q'R'P'. 
Examining the diagonal elements of R' we can determine its rank r'; the columns of 
R' correspond to a permuted basis of Jk H (A/JjjLi) and the generators corresponding 
to the first r' columns are redundant and can be discarded. 

Exact Gaussian elimination, SVD and QRP-algorithm applied to an m X n 
matrix all have the same complexity 0(mn min(?7j, n)). In the approximate algo- 
rithm the computation of the null-space using SVD has a complexity 0(nh 3 ) and is 
followed by a QRP-algorithm which has a complexity 0{n 3 h 3 ). Thus the complex- 
ity to compute a set of generators or a border basis of I{J\, . . ■ ,J S ) is 0(sn 3 h 3 ), 
while the complexity to compute a minimal set of generators of I{J\, ■ ■ ■ , J s ) by the 
algorithm MinimalBasis is 0(sn 6 h 3 ). 

Alternatively, after computing a border basis, we can give a numerical algorithm 
based on Proposition 12 .141 which would reduce the complexity to 0(sn 4 h 3 ) with a 
slight loss of numerical precision. 

Example 4.1. We implemented our algorithm in Octave. Here is the result 
we obtained when we tested it on the following parametric sextic space curve C 
taken from [JWGJ: 

x = 3sH 2 - 9s 3 t 3 - 3s 2 t 4 + 12st 5 + 6t 6 
y = -3s 6 + I8s 5 t - 27sH 2 - I2s 3 t 3 + 33s 2 i 4 + 6st 5 - 6t 6 
z= s 6 - 6s 5 t + 13sH 2 - 16s 3 t 3 + 9s 2 t 4 + Ust 5 - 6t 6 
w = -2sH 2 + 8s 3 t 3 - Us 2 t 4 + 20st 5 - 6t 6 . 
By Theorem 15.11 the ideal of this curve of degree 6 in P 3 (C) can be generated 
by polynomials of degree at most 5. 

We chose 31 > 6 • 5 points using roots of unity of the following form: 

s = l; t = exp(2Tdk/31) k G {1,...,31}. 

Running the algorithm BorderBasisWithComplement, we obtained no polyno- 
mials of degree 1 (showing that the ideal is not contained in any hyperplane), no 
polynomials of degree 2, 4 polynomials of degree 3, 11 of degree 4 and 22 of degree 
5, yielding a set of generators for the ideal in .08 seconds. Among the 20 monomials 
of degree 3 the QRP-algorithm chose z 2 x, yxw, y 2 w, z 2 w as border monomials and 
the remaining 16 as generators of a complement. 

We obtained a minimal basis consisting of only the 4 polynomials of degree 3 
in an additional time of .02 seconds: 

/i = z 2 x + 0.06666666666z 2 ?/ + 0.06805555555 zy 2 - 0.0361111111 zyz 

-0.2833333333 zyw - 0.55 zx 2 - 1. 066666667 zxw + 0.01527777778y 3 
-0.09166666666y 2 x + 0.3055555555 yx 2 + 0.1833333334a; 2 w, 

h = yxw + 0.2 z 2 y + 0.1416666667zy 2 - 0.4833333333 zyx - 0.1 zyw 
-0.9 zx 2 - 0.2000000001 zxw + 0.025 y 3 - 0.15 y 2 x + 0.5 yx 2 
+0.3000000001 x 2 w, 
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h = y 2 w - 0.8 z 2 y - 0.3166666667zy 2 - 0.5666666667 zyx + OAzyw 
+0.6000000001 zx 2 + 0.8000000002 zxw - 0.01666666666y 3 
+0.09999999999 y 2 x - 0.3333333333 yx 2 - 0.2000000002 x 2 w, 

f 4 = z 2 w - 0.6666666667z 3 - 0.162962963 z 2 y + 0.06049382717zy 2 
-0.03209876545 zyx + 0.9703703704zyw - 0.4888888888 zx 2 
+0.3851851853 zxw - 0. 1666666667 zw 2 + 0.01358024691 y 3 
-0.08148148149 y 2 x + 0.2716049383 yx 2 - 0.9444444445 yw 2 
+0. 1629629629 x 2 w - 0.2222222224xw 2 

Using continued fractions, we then attempted to convert the coefficients from 
floating point numbers to rational numbers, obtaining the following polynomials: 

fi - 7 2 x + J- z 2 7, 4- -*?- 77; 2 i il „3 _ ii 71IT _ JJ_ ,.2 T _ 11 7T 2 
j I — ii a t 1 - * y T 72Q ^.y T 720 y 360 120 » 20 

11 9 17 16 119 

+ 36 V X ~ 60 Z V W ~ 15 ZXW + 60 ^ W ' 

/ 2 = yxw + i z 2 y + ^ zy 2 + i y 3 - §§ zyx - ^ y 2 x - ^ zx 2 + i yx 2 

— zyw — i zxw + ^ x 2 w, 

/ 3 = y 2 w _ 4 z 2 y _ 19 zy 2 _ i y 3 _ II zyx + 1_ y 2 x + 3 ^2 _ 1 ^2 

+ | zyu> + I zxw — I x 2 w, 

^ = " I " 3 22 ^2 Z \ + 1 ^ 22/2 5 t ^ ^ " 2 F f ^ "l 1 ^ f X 17 2 

-|g zx + ^ yx 2 + zyw + f^ zxw + x^w - i zw 2 - j± yw 2 

2 9 ' ' 

— I HO . 

The floating point coefficients were sufficiently accurate to recover the exact 
rational coefficients and the previous polynomials generate the exact ideal of the 
curve over Q. 



5. Degree bounds for ideal generators 

Let C be an irreducible projective curve (seen as a set of points in P"(K)) and 
let / = 1(C) C V — K[xo, . . . , x n ] be the prime homogeneous ideal consisting of all 
the polynomials vanishing on C . In Section [3] we saw that the computation of / 
can be reduced to the computation of the ideal of sufficiently many points on the 
curve (see Proposition 13. ip . In order to be sure that one has enough points, it is 
necessary to bound the degrees of generators of the ideal /. Such a bound can be 
obtained from the regularity of the curve. 

Recall that if — > . . . — > F\ — > F — > J — s- is a graded free resolution of 
a polynomial ideal J , we say that J is fc-regular if each Fj can be generated by 
polynomials of degree < fc + j. Then we call regularity of J the integer reg(J) = 
min{fc | J is fc-regular}. 

If reg(C) = reg(X(C)) — m, then (see for instance [E]) I can be generated 
by homogeneous polynomials of degree at most m; moreover the Hilbcrt function 
HFi(s) and the Hilbert polynomial HPi(s) of V /I take the same value for all 
integers s > m. 

The next result gives information about the regularity of the curve as a function 
of the curve degree and of the dimension of the embedding projective space: 

Theorem 5.1 (Gruson-Lazarsfeld-Peskine |GLPj ). Let V C P n (K) he a re- 
duced and irreducible curve of degree d not contained in any hyperplane, with IK 
algebraically closed and n > 3. Then reg(T>) < d — n + 2. 
Moreover, if V has genus g > 1 then reg(V) < d — n + 1. 
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The previous result and many of the degree bounds apply to curves not con- 
tained in any hyperplane, i.e. non-degenerate, which happens if and only if I\ = 
1(C) i = (0). On the other hand a degenerate curve in P n (K) is a non-degenerate 
curve in the projective subspace V(I±) having dimension n — dimii. So the bounds 
of this section apply to degenerate curves if we replace n by n — dim Ji . 

Remark 5.2. A parametrization of a curve C C P"(K) can be regarded 
as a machine delivering as many points on the curve as needed. Moreover if 
%l>: P 1 (K) — > ¥ n (K) is a rational map whose image is the curve C and ip([to,ti]) = 
[Fo{to, ti), . . . , F n (to, ti)] with Fi(to, ti) homogeneous polynomials of degree s, then 
degC < s. Using this upper bound for the curve degree, via Theorem 15.11 we ob- 
tain an upper bound for the degrees of generators of C. This allows us to use our 
procedure to compute generators of the ideal / of the curve and in particular to 
compute an implicit representation of C starting from a parametric one. □ 

Sharper bounds for the degrees of generators can be obtained for certain curves 
obtained as the image of the embedding given by a complete linear series. A first 
result in this direction concerns canonical curves (see for instance |SD2j ): 

Theorem 5.3 (Petri [P]). The ideal of a non-hyperelliptic canonical curve of 
genus g > 4 can be generated by polynomials of degree 2 and of degree 3. 

Proposition 5.4. Let C C P"(K) be a non-degenerate irreducible curve of 
geometric genus g > 4 and of degree d = 2g — 2, with K algebraically closed and 
n = g — 1 . Then the ideal / = 1(C) can be generated by polynomials of degree 2 
and of degree 3. 

Proof. By Theorem VI. 6. 10 in |W] any non-degenerate curve of genus g and 
degree 2g — 2 embedded in P 9_1 (K) is a smooth non-hyperelliptic canonical curve. 
Then the conclusion follows from Theorem 15.31 □ 

The following result of Saint-Donat gives bounds for the degrees of generators 
of curves of genus g which are the image of an embedding given by a complete 
linear series of sufficiently high degree: 

Theorem 5.5. ^ [SDl] ) Let C C P n (K) be an irreducible nonsingular curve of 
genus g embedded by a complete linear series of degree d. 

(1) // d > 2g + 1, then the ideal I of C can be generated by polynomials of 
degree 2 and of degree 3. 

(2) If d > 2g + 2, then I can be generated by polynomials of degree 2. 

The following result shows that we can remove the condition of being embedded 
by a complete linear series: 

Proposition 5.6. Let C C P"(K) be a non-degenerate irreducible curve of 
geometric genus g and degree d, with IK algebraically closed and n = d — g. Then 

(1) If d > 2g + 1, then the ideal / = 1(C) can be generated by polynomials 
of degree 2 and of degree 3. 

(2) If d > 2g + 2, then / = 1(C) can be generated by polynomials of degree 
2. 

Proof. Recall that any curve C of degree d and genus g in P n (K) can be seen 
as the image of the map given by a linear series contained in C(D) for some divisor 
D of degree d on a non-singular model C of C . 
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For any divisor D on C, denote by 1(D) the dimension of the complete linear 
series C(D). We also denote by K a canonical divisor. 

If cZ > 2g + 1, then l(K — D) = and hence by Riemann-Roch Theorem 
1(D) = d — g + 1; moreover (see [F]) the map induced by D is an embedding. 
Thus we can see C as the image of an embedding of C in P d "f(K). So, since C is 
non-degenerate, C is embedded by a complete linear series if and only if n = d — g. 

Moreover, if d > 2g+ 1 and n = d—g, the curve C is the image of a non-singular 
curve through a map which is an embedding, therefore C is non-singular. Hence 
the hypotheses of Theorem 15.51 are fulfilled and it implies our result. □ 

Example 1. Assume that C has genus g > 3 and is embedded by the bicanon- 
ical map. Since deg(2K) = 4g — 4 > 2g + 2, by Theorem 15. 51 / can be generated by 
quadratic polynomials. 

Example 2. If C has genus g — 2 and is embedded by the tricanonical 
map (which gives an embedding because deg(3K) — 6g — 6 — 6 > 2g + 1), since 
deg(3AT) > 2g + 2 again by Theorem 15.51 / can be generated by polynomials of 
degree 2. 

Other results giving bounds for the degrees of generators of the curve ideal in 
different situations are available in the literature; the following one (see [A]) gives 
results in the hyperelliptic case: 

Proposition 5.7 (Akahori). Let C c P"(K) be an irreducible non-degenerate 
and non-singular hyperelliptic curve of genus g and degree d. 

(1) If d = 2g, then the ideal I oiC can be generated by polynomials of degrees 
2, 3 and 4. 

(2) If d = 2g — 1, then / can be generated by polynomials of degrees 2, 3, 4 
and 5. 

When we need to compute the null-space of the evaluation matrix Mr for 
approximate points (see Section [4]), we can try to determine its rank by inspecting 
its singular value spectrum. Since this is not guaranteed to work, it is useful to 
predict what the rank should be. Since the rank of M-ji equals |A4| = dimP^ — 
dim Jfe, the following proposition gives a situation where we can predict this value. 

Recall that if go, . . . , g n is a basis of the complete linear series C(D) of dimension 
1(D) =rt + l, then Ik = Kerifk where ipk- Vk — > C(kD) defined by tpk(%i) = 9i- 

Proposition 5.8. If tp k is surjective and k ■ deg(-D) > 2g — 1, then 
dim"Pfc — dim/fe = k ■ deg(D) — g + 1. 

Proof. If ipk is surjective, then dim Ik = dimKer^fc = dim 7-^ — dim C(kD). 
Since deg(kD) = k ■ deg(D) > 2g — 1, then i(kD) — 0. Hence by Riemann-Roch 
Theorem we get l(kD) = deg(kD) — g + 1 = k ■ deg(D) —g + 1 and the conclusion 
follows. □ 

Examples. 1. Assume that C is non- hyperelliptic of genus g > 4 and take 
the canonical divisor D — K; in particular deg(£>) — 2g — 2 and 1(D) = g (i.e. in 
this case n = g — 1). By Theorem 15.31 we already know that / can be generated 
by polynomials of degree 2 and of degree 3, and so it suffices to know dim/2 and 
dim ^3 . 

Then by Noether's Theorem the map ifk- Vk — > C(kD) is surjective for all k. 
Furthermore for all k > 2 we have that k ■ deg(D) > 2g — 1, hence by Proposition 
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15.81 we have dim Vk 



dim/fe = k(2g — 2) — g + 1. In particular 
dimP 2 - dim J 2 = 2{2g -2)-g+l = 3g 
dimP 3 - dim/ 3 = 3{2g -2)-g + l = 5g 



5. 



3 



2. If we choose a divisor D such that deg(D) > 2g + 1, then (/3j. : — > C(kD) 
is surjective (see |Mj); moreover for all fc > 1 we have that k ■ deg(D) > 2g — 1. 
Then we can compute dim 1^ by Proposition 15.81 

In particular if D = 2K and g > 3, then we know that I can be generated 
by quadratic polynomials. Since i(2K) = 0, by Riemann-Roch Theorem l(2K) = 
3g — 3, i.e. n — 3g — 4. Then 



If g = 2 and we choose D = 3K, we know that / can be generated by quadratic 
polynomials. Since i(3K) — 0, by Riemann-Roch Theorem l(3K) = 5g — 5 = 5, i.e. 
n = 4. Then 



Remark 5.9. When the ideal can be generated in degree 2 and we know |A/2| 
as above, then the algorithm to compute a minimal set of generators can be consid- 
erably simplified. If M-jz is the point evaluation matrix for all monomials of degree 
2, since we know that rk M-jz = IW2I, a single application of the SVD-algorithm 
with this imposed rank computes a basis for the null-space of Mr which directly 
yields a minimal set of generators of the ideal. 



We wish to thank Mika Seppala, for proposing this problem to us and a referee 
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